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ABSTRACT 

This is the third paper on the improvements of systematic errors in our weak 
lensing analysis using an elliptical weight function, called E-HOLICs. In the 
previous papers we have succeeded in avoiding error which depends on ellipticity 
of background image. In this paper, we investigate the systematic error which 
depends on signal to noise ratio of background image. We find that the origin 
of the error is the random count noise which comes from Poisson noise of sky 
counts. Random count noise makes additional moments and centroid shift error, 
and those 1st orders are canceled in averaging, but 2nd orders are not canceled. 
We derived the equations which corrects these effects in measuring moments and 
ellipticity of the image and test their validity using simulation image. We find 
that the systematic error becomes less than 1% in the measured ellipticity for 
objects with S/N > 3. 



Introduction 



The importance of the weak lensing analysis is now widely recognized because it has a 
potential to provides us a direct and unbiased information on the mass distribution for lens 
objects. The weak lensing analysis measures shapes(called ellipticity which has two com- 
ponents interpreted as direction and magnitude) of many background images (galaxies) and 
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then averaged over an appropriate number of images to get rid of intrinsic random ellipticity 
of images and to withdraw the elhpticity due to gravitational tidal effect (shear) of the lensing 
object. The shear carries the information of mass structure of the lensing object. Thus an 
accurate shape measurement of the background images is critically important to accurately 
measure the mass distribution. So far weak lensing is very successful for cluster lensing (ellip- 
ticity due to shear is of the order of 5%) and provides us a rich information of mass structures 
of clusters and of our understanding structure formation in the universe. Recently the cos- 
mic shear, i.e. the weak lensing due to large scale structure(LSS) attracted much attention 
because of it's ability to study the nature of dark energy which is supposed to be the source 
of the accelerated expansion of the universe. In fact several projects for the cosmic shear 
measurement are proposed and some of them is almost ready to start the observation (Hyper 
Suprime-Cam http:/ /www.naoj.org/Projects/HSC/HSCProject.html, Dark Energy Survey 
http:/ /www. darkenergysurvey.org/, Euclid http:/ /sci.esa.int/euclid and so on). However the 
signal of cosmic shear is very weak (of the order of 1%) compared with cluster lensing and 
thus needs special treatment. Namely we needs to develop very accurate shape measurement 
scheme which avoids various systematic errors. For example the measured gravitational shear 
depends on the ellipticity and signal to noise of background image. Usually such dependence 
becomes small by averaging many of the images, but it is critically important to realize that 
these dependence somehow correlated with the redshift distribution of image which is also 
important to have an accurate measurement of the shear. Thus we cannot make a simple 
averaging over the images without having a method free from such systematic errors. The 
required accuracy for the measurement of ellipticity is less than 1% in order to have an useful 
information of dark energy. 

There have been many studies in this direction and various measurement schemes are 
proposed(Kaiser et al 1995, Bernstein &; Jarvis 2002; Refregier 2003; Kuijken et al. 2006; 
Miller et al. 2007; Kitching et al. 2008; Melchior 2011). The accuracy of these methods are 
tested using the simulation data provided by STEPl(Heymans et al 2006), STEP2(Massey 
et al 2007), GREAT08 (Bridle et al 2010) and GREATIO (Kitching et al 2012). Although 
much progress is reported, none of the method achieved the required accuracy and are free 
from various systematic errors. 

We have also developed a new scheme based on KSB method(Kaiser et al 1995) using 
an elliptical window function(we called E-HOLICs) to measure the background image as 
accurately as possible (Okura and Futamase 2011, Part I paper). It is shown in our Part II 
papers that the E-HOLICs can improve the systematic error which depends on ellipticity. 
In this paper we study the systematic error which depends on signal to noise ratio(SN). 
There are some studies about this systematic error (Hirata et al 2004, Kacprzak et al 2012, 
Refregier et al 2012, Okura and Futamase 2012 Part II paper, Melchior and Viola 2012). 
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these results show this systematic error comes from random count noise(RCN). Because, 1st 
order effects from RCN arc canceled by averaging, but 2nd order effects are not canceled. 
We calculate the 2nd order effects to obtain the correction formulas in the measurement 
of moments and ellipticity for Gaussian weighted images in KSB method and E-HOLlCs 
method (i.e. without PSF correction). Using the simulation data GREAT08, we find that 
the derived formula correct the SN dependent bias within 1% for images with SN > 3. 

The paper organized as follows. In section 2, we explain and define our notations and 
some of the definitions. In section 3, we calculate the 2nd order effects of RCN and obtain 
general formulas to correct the effects. We test the formula in the case of KSB method with 
Gaussian weight function. The correction formula in the case of E-HOLICs is presented 
in section 4, and tested it using GREAT 08 simulation and find that the systematic error 
becomes less than 1% in the measured eUipticity for objects with S/N > 3. In section 5, we 
summaries our results. 



2. Basis and Definitions 

In this section, we present notations and definitions we use in E-HOLICs method. Some 
of them were defined in part2, but we add the effect of "random count noise" (hereafter RCN) 
and " centroid shift error" (hereafter CSE) . 



2.1. Random Count Noise 

First, we write the observed brightness distribution of object as " I°^^[d)" , which is the 
sum of object "/''^J(^)" and RCN "/^^^(e)", so 

where, "0" is position angle in complex coordinate whose origin is at the centroid of object 

e = el = ei + i02, (2) 

and the products of the positions are notated as 

N+M N-M 

oz - (01)^ (3) 

N means order and M means spin-number. 
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We assume that RCN is Poisson noise of sky counts and also assume that all pixels have 
same root-mean-square(RMS) of RCN "ctjicn" ■ We don't consider Poisson noise of objects 
itself in this paper. If an object has a photon count Nobj, then the Poisson noise is of the 
order of and thus the order of errors reduces by 1/ ^jNobj- Therefore if the object 

is bright enough to be able to neglect sky noise, we can also neglect own Poisson noise. On 
the other hand, if an object is faint, its Poisson noise is much smaller than the Poisson noise 
from sky Nsky{Nobj ^ Ngky), so we can neglect it. However it needs another consideration 
for the situation with Nobj ~ Ng^y which will be discussed in other paper. 



2.2. Notations 

In measuring moments of image by E-HOLICs method, we use an elliptical Gaussian 
weight function with ellipticity " dw = ^wi + iSw2" for measuring the complex moments, 
and we define this weight function as 



W{0,Sw) = e ^ , (4) 
where cx^ is a size parameter of weight function. 

The complex moments and HOLICs of arbitrary brightness distribution without centroid 
shift error(CSE) are defined as 

zZii,Sw) ^ j d^eel,i{e)w{e,5w) (5) 

UN/T ryO X \ — ^ui^l^w) (n\ 

In this paper, we define the origin of the coordinate at the centroid of I°^^{0), therefore 

zl{r'\Sw) = 0. (7) 

However, RCN causes CSE, so we cannot obtain eq.(l7l) in real analysis. We notate CSE due 
to RCN as A6 = A6l. then the complex moments with RCN and CSE that we measure in 
real analysis are defined as 

and HOLICs are measured as 



yN ( jobs s \ 

njNtTobs yO s \ — ^ i^w) /^.n 



The detail of this CSE is expressed in section 3. 
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2.3. WSN 

Here, we define weighted signal to noise ratio "WSN" with elhptical weight function as 
ctrcnJ J dWie^Sw) (^RCN\/S^ (TrcnVS^ ' 



where Sw is an integral of weight function or weighted area 



5w = \Sw\- (12) 
WSN appears frequently in the following calculations, so we use WSN instead of SN. 

We measure SN and WSN of back ground objects detected from Abell 1689 real data 



taken by Subaru suprime-cam, and we use only objects having > 7 by IMCAT(http: //www.ifa.hawaii.edu/ 
detection. And we define a signal to noise ratio "SNj," defined from u as 

SiV„.^. (13) 

The plots of SNu and WSN are shown in figH] and we can see the following relation. 

WSN ^ 3SN^. (14) 

Figl2] shows the count distributions of SN,y and WSN. 



2.4. Averaging 

In weak lensing analysis, usually "averaging" means the averaging over a parameter of 
several different objects, but in this paper "averaging" means averaging over a parameter 
of same object but different RCN. It means that we observe the same object many times. 
Therefore the difference between the value measured without RCN and the averaged value 
with RCN is the systematic error. 

Because we use different RCN, the averaged complex moments of them vanish, 

Z^il^cN^Sw) = 0. (15) 

So averaged values of complex moments of RCN are 0, but squares of the moments are not 
0, because RCN has self-correlation. 
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Fig. 1. — Plots of SN,, .v.s. WSN. We used background objects of Abell 1689 real data. We 
reject objects which have SN lower than 2. We can see some objects have correlation about 
WSN = 3SK and almost objects have WSN ^ 3SN^. 
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Fig. 2. — Distributions of counts of SN and WSN. We used background objects of Abell 
1689 real data. We reject objects which have SN lower than 2. 
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Standard deviation of RCN at arbitrary position Oa is obtained as 
1 ^ 

i 

where "i" means ith set of RCN and 6d{0) is Dirac delta Function.. Let G^^ be defined as 
a square of complex moments of I^^^ {6) as 



yN(TRCN s \yO{TRCN s \ 



^RCN 



and average value of G^j is obtained as 

1 



^RCN 

1 

RCN 

1 

° RCN 



d^eeZwie, dw) [ d^9'e'^w{e\ 5w)5d {o - e') 



= j d^eeZX%{w{e,dw)f. (is) 

Because we use elliptical Gaussian for weight function W{6, 5w)i G^m(^vk) can be calculated 
analytically. The detailed values of GZ{5\y) can be seen in Appendix lA.li Then we can 
obtain the averaged value of the square of the complex moments as the product of c'^cn ^^cl 
Gm{Sw)- 

Here we show explicitly the calculation for a\{I^'-^^ , ^w)\+ which is the standard devi- 
ation of each components of ZKl^*-^^ , Sw),- 



' '^RCN 



I d^e{e,w{e,dw)f + iJalcNj dW{e2W{0,5w)y 



\ 2^ ^ + 1 



Glidw)-Re Gl{6w 



( 

<^RCN 

\ 

figj3] shows simulation results of eq.( IT9|) with normalization by ^,/Swo^j2 and aRCN = 1- 
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Fig. 3. — Test of standard deviations of Zl{I^'-'^ ,5w) with normalization by ^/Swo^^ 
and aRCN = 1- Horizontal axis means ellipticity of weight function Swi (ySw2 = 0) and 
vertical axis means a standard deviation of the complex moments. Cross(square) points mean 
real(imaginal) part of al and line(dash) means real(imaginal) part of analytical estimation 
eq.din]) 
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3. Centroid Shift Error and Complex Moments with Random Count Noise 

In this section, we present calculations of centroid shift error(CSE) and the complex mo- 
ments with random count noise(RCN) in detail, where we assume the brightness distribution 
of object as an elliptical Gaussian image and adopt KSB method. 



3.1. Centroid Shift Error with Random Count Noise 

We present calculations about CSE. In this paper, we define true centroid as the origin 
of the complex coordinate, so = is true centroid, however centroid of we measure 

is different from the origin and this difference is CSE " A0" . 

We measure the centroid as a position which vanishes the dipole moment of I°'^'^{d), so 
the dipole moments of I"''^{0) is 

zl{r'',dw) = zl{r'^,dw) + zl{i''^'',Sw)^o, (20) 

and by expanding with AO and neglecting higher order, we obtain 

zl{r'% dw) = J d^e {6 - Ae)l r'^{e)w{e - ao, dw) + zKi''^'', Sw) 

{r'\dw) + zl{i^'''',dw) 



W{Ad,6w) 



^ _ W{Ad,8w)Z'o{I'"',Sw) ^^20^ + cot Ar) + Zl{I^^\ Sw) = 0, 



(21) 



where we use the elliptical Gaussian form for the weight function to write the weight function 
in the following form 



/ (A9\* -d*wA9l)el + (A9l-dwA9l*)9l*\ , 



A9§-Re[g;^A9|] 

W{Ad,6w) = e ^ , (23) 
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C coefficients are defined as 



nN 



pN 



'W 



a 



w 



- {{N - M) n^-' + {N + M-2) SwUZzl) 

I {{N + M) n^-' + {N-M-2) s*^nZ-^l) - 2- 



1 - 

^ 2 

'W 



•w 



w 



(njN+2 c* njN+2\ 



{r'^,z^,m 

(26) 
(27) 



These have non dimension and have spin-M, X is an integer. Here we neglect odd order of 
the complex moments of r''^{d) (i.e. Zl{r''\5w) = Zl^{r''\6w) = 0). Finally, we obtain 
AO as 



2 


C2l_Zl - COl_Zl* 




RCN s \ 


Z^,{lo^^,Sw) \C2^o+ 


2 - |COt 


2 



Zl{I'^'^,6w) 



C2l_Z\-C^\_Z\* (/^^^,5i^), (28) 



where 



r2° = 



l/^oO 12 



|C0t|2 



(29) 
(30) 



Therefore, Ad comes from the dipole moments of RCN. 

The averaged value of Ad is 0, because averaged value of Z\{1^'^^ is also 0, 
therefore there is no SCE in averaged sense. But the averaged value of (A^)^ is not 
and obtained as follows. 



w{Ae,8w)Ae = 



W{AQ,bw)Ael 



W\Ae,5w)AQ 



1 


2 


[\C2l\' 


WSN"^ 1 




1 


2 




WSN^ 1 





2Re 



C2',_C0tdw 



We define ellipticity of distribution of CSE as 



Aei W{A0, dw)Aet, 



Aei w{Ae, 5w)Aei 



Sc = \S, 



c\ 



(31) 

(32) 

(33) 



(34) 
(35) 
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It the ellipticity due to CSE distribution. 

If the object has an elliptical Gaussian image I^^{6, dj) 



/^^(6>,5,) = (36) 

where A is an arbitrary amplitude (hereafter EGI) and we use KSB method with Gaussian 
weight (hereafter KSBGW). We can obtain eq.f p2|) and eq.fl33l) analytically as 

A9l^W^{Ae,S^)A9l = + (37) 

3 



WSN^ 








WSN2 






3X2 
2" I 




5X2 
4" I 



Ael^W\Ae,6w)Ael = ^^{l + -6j] (38) 



= Si^-^^[^ + \sjjdr, (39) 

where the detail value of C coefficients in this situation are shown in Appendix IA.2I 

FigJH shows the comparison between the result of simulation and the predicted result 
using eg. (1371) and eq.(l38l) with Sj = (0.5, 0) and o"^ = 200. We can see these equations are 
very good approximation except in low WSN. 



3.2. Complex Moments with Random Count Noise 

In this section, we consider the systematic error in measuring the complex moments due 
to RCN. 

By expanding eq.([HD, Z^^ can be expressed as a function of Z^, so 
ZZi^'^ Sw) = J d'e {9 - A^);, F'^{d)W {0 - AO, 6w) 

/(Ae*-«^Ae)e+(Ae-6T4/Ae*)e* 
(fe {9 - A9)^ r''{0)W (6>, Sw) W {A0, dw) e ^ 

= / d'9 (9 - A9)^, (r'^iO) + I^^^iO)) W (6, Sw) (l - - M^w^^l] 



, ^ , {AO* - S*^A0) 6 + (AO- SwAO*) 0* , 1 ,2 
X I 1 + -2 + — ^ [AO - S^AO) 92 

+-V (A^ - SwA0*f 9f + \{A0~ SwA0*)l 9^ + ...]. (40) 
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WSN 



Fig. 4.— Plots are W^AO, 6w)Ael^/aw in the situation of EGI and KSBGW. Hor- 
izontal axis means WSN and vertical axis means W'^{A6,Sw)'A9l^/a^. Square(circle) 
means W'^{A6,dw)A9Q/a^{W'^{A0,Sw)A92/(j'^) measured from simulation data and 
solid(daslied) line means analytical prediction i.e. eq. fl32|) (eq.f l33|) ). 
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Here, we define the effects of RCN in measuring tfie complex moments as 

5w) = z^{r'\ 5w) [nUi"''. z^, ^w) + ^n'^Mir^idw) + ^hI(2){6w)\ , (4i) 

wliere X of ^l-L^f^-^-^ means from Xtli order effect of I^'~''^[6). Because effect from 1st order 
of /«^^(6») is 0, so 



A^M(i)(^^) = (42) 



AnZ^,^{Sw) ^ (43) 



A'HZ^,^iSw)AZ^^^)iSw) ^ 0. (44) 
1st and 2nd order of effects for complex moments from I^'-'^[0) can be written as 



where 



''2 _ (^2 ^*-'M-2+ *-''JM-2+ J ^ ^2 ^^"1^ _ ^2 ^ *-'M+2- *-''Jm+2- 

+-"2- fc2o^G'j^+^ - C02lGj^+2 ~ ^iV^2o_Gjlf^2^ + ^^C02_G'm"^^) ) , (46) 



eg = A^o' - '5;^A^2' (47) 
= A9l-dwAel (48) 



So, we can obtain systematic error ratio " A?^jJ^^2)(^i^)" by C coefficients (i.e. combinations 
of complex moments) . 

Here, we define systematic error ratio (hereafter SER) of complex moments as 
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and 2nd order of systematic error ratio(hereafter 2ndSER) of complex moments as 



If 2nd order effect of RCN is dominant, these SERs are almost same(SER?«2ndSER). So, a 
correction formula for HOLICs which corrects systematic error is defined as 



^ HUi'"\z^,Sw)-^Km^Sw). (51) 

In the situation of EGI and KSBGW, we can obtain the average effects of RCN as 

1 2 + Q, 



^^0(2) (0) - WSN^ 4 ^ o(/^^^o°,0) (52) 

"^"I u2f tEG 



^^0(2) (0) ^ -^SF 8 ^^^^ 

A^2^2)(0) ^ -^§^0 + '^')^^'^^'''''^°'°^ ^^^^ 



1 /I 3 



^^0(2) (0) - 2 + 4^^ro(^ 



1 7 



^^2(2)(o) - "WsFi^ + s^V^^^^ ^^^^ 



1 /I 9 



^^4(2)(0) - -WsFl^2 + 4^^j^'^^ ^^^^ 

Figl5] to FigJH] show simulated results of SER and 2ndSER calculated from these equations. 
Fig 15] shows about monopole moments and eq. (l52|) . Fig 16] shows about spin-0 quadrupole 
moments and eq. (l53]) . FiglT] shows about spin-2 quadrupole moments and eq. (l5^ . FiglH] 
shows about 16pole moments and eq. (15Bl) . 



3.3. Ellipticity with Random Count Noise 

Now we show the systematic error in measuring of ellipticity due to RCN which is used 
in the weak lensing analysis. The observed ellipticity and object ellipticity are defined as 

dobs = n2[I , Zq, dvy) = — — — — (58) 



= 'Hl{r'^,Zl6w)= ^l)r.J,i , (59) 
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1 10 

WSN 

Fig. 5. — SER due to RCN for monopole moments measured in the situation of EGI and 
KSBGW. Horizontal axis means WSN and vertical axis means SER. Plots are simulated 
results of SER of AT/q and line is 2ndSER calculated from eg. (1521) . 
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Fig. 6.— SER due to RCN for Hi in the situation of EGI and KSBGW. Horizontal axis 
means WSN and vertical axis means SER. Plots are simulated results of SER of A'Hq and 
line is 2ndSER calculated from eq.( l53|) . 
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1 10 

WSN 

Fig. 7. — SER due to RCN for measured in the situation of EGI and KSBGW. Horizontal 
axis means WSN and vertical axis means SER. Plots are simulated results of SER of A'H2 
and line is 2ndSER calculated from eg. (15^. 
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WSN 

Fig. 8.— SER due to RCN for nj^ in the situation of EGI and KSBGW. Horizontal 
axis means WSN and vertical axis means SER. Square (circle, triangle) plots are simulated 
results of A'Ho(A'H|, AT-Lf) and solid (long dashed, dashed) line is 2ndSER calculated from 
eq.(l55D(eq.(l56D,eq.([S 
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and we consider effect from RCN until 2nd order as 

^obs ~ ^obj 

+ A<5(i) + A(5(2) (60) 

wfiere number of sub scripts mean order of RCN effect. Here, we define systematic error 
ratio (liereafter SER) of elliticity as 



SER = ^"'^ (61) 

and 2nd order of systematic error ratio (liereafter 2ndSER) of ellipticity as 



2ndSER = (62) 

Oobj 

If 2nd order effect of RCN is dominant, these SERs are almost same(SER?^2ndSER) So, a 
correction formula for ellipticity is defined as 

^obj (corrected) ^ ^obs _ j^gg^j 

For such combinations of complex moments, we cannot neglect the 1st order effects 
because there are combinations of the 1st order effects. For example, RCN effects for the 
observed ellipticity can be calculated as 



'obs 



dobj + ^'^2{l)i^w) + A'H2(2)(<5vi/) — <5ofejA'Ho(i)(<5vi/) — Sobj^'Hl(2){^w) 

{Hl{P^^,Zldw))\Gi-5wGl) 



w 



(64) 



so 



A(5(i) = Anl^,^{dw) - SobA'^mi^w) (65) 
Ad(2) = /XH2{2)[Sw) ~ OobjAnQ^2)[^w) ^^^^^^^2 (66) 



In the situation of EGI and KSBGW, the averaged value of eq. flMj) can be calculated 
analytically as follows. 



^obs 



WSN" V2 

obj /I 9r2 



AS(2) = + (68) 

^ ' WSN^ V 2 8 ^ ' ^ ' 
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Fig|9]is the comparison between the resuh of simulation and prediction by eq. fl68|) with di = 
(0.5, 0) and in this test we use the same size of image for weight function, so dobj = (0.25, 0). 
We can see eq.f pK]) gives a good agreement except for the sources with low WSN, and the 
ellipticity of object having WSN = 5 is measured with about 4% underestimation. 



4. E-HOLICs method with Random Count Noise 

E-HOLICs method uses the ellipticity of objects for the ellipticity of the weight function. 
However we cannot measure true ellipticity due to random count noise(RCN). Therefore we 
must consider RCN effect for the ellipticity of weight function. 



4.1. E-HOLICs method with Random Count Noise and true ellipticity for 

weight function 

In this section, before treating the realistic situation, we consider the ideal situation 
where we use true ellipticity for the weight function. Thus the ellipticity is written as 
follows: 

S.. ^ (69) 
Zi{lobj,Sw) ZiilobJ^Sj 



^Obj - ry2tTohi X \ - ry2tTohi X \ - - ^1 (70) 



Calculations in this section are useful for consideration for the realistic situation which is 
explained in the next section. C coefficients of elliptical Gaussian image (EGI) are calculated 
analytically and shown in Appendix IA.3I 



Centroid shift error (CSE) for the image with arbitrary distribution is given by eq. (!28|) . 
and average effects are given by eq. (l32|) and eq. (l33|) . 

Ae^WiA0,dw)A0 = I \C21_ZI - CQlJl*] (J^^^,^,), (71) 

Zo[I"''^,Oi) L J 



Ael^W^{Ae,Sw)Ael = 1 1\ + \COU' - 2Re [C2l_C0tS^]) (72) 

Ael^w^{Ae,dw)Ael = (^(c2°_)'^, -2C2°_co°_ + (cot) '^1^73) 
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Fig. 9. — SER due to RCN for observed ellipticity in the situation of EGI and KSBGW. 
Horizontal axis means WSN and vertical axis means SER A5(2)/<5. Plots are simulated 
results and line is 2ndSER calculated from eg. (168!) . 
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The average CSE effects of EGI can be obtained as 

M^W\^e^6w)Ml = (74) 

Sc = Si. (76) 

FigJTOl shows the comparison between the results of simulation and the prediction be eq.( !74|) 
and eq.([75D with 6i = (0.5,0). 

From eq.dMD AH^^(2) 

are obtained as 
/iV + M _ A , _ 

\^^^ _ ^2^^1-2+ ^^M-2+J^^2 y^^^ _ ^2-s^^M+2- ^^M+2- 

+^ (C2°^Gjr/^ - cor - S*jC2lGZl + ^^^COLGj^/^) ) , (77) 

The average effects for the complex moments of EGI can be obtained as 

'm^) - ^^^o(/^^,^o°,5/) (78) 
^ (79) 
AH2^(2)(5/) ^ (80) 

^ ^^^^ 
AH^(,)(5/) - -^^^^(/^^,^o^'5.) (82) 

FigHH to figll3] are the comparison between numerical results of SER and theoretical pre- 
dictions of 2ndSER calculated from eq. flTSj) to eq. fl83|) . respectively. All cases we assume 
(5, = (0.5,0). 
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Fig. 10. — Plots are A6'^/(T^ with elliptical Gaussian image <5/ = (0.5,0) and E- 
HOLICs method dw = ^i- Horizontal axis means WSN and vertical axis means AOj^/a^. 
Square(circle) means A6'Q/cr^(A6'|/cr^) measured from simulation data and sohd(dashed) 
line means analytical prediction i.e. eq.f l7^ (eq. fl75|) ). 
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Fig. 11. — SER due to RCN for monopole moments measured with EGI and E-HOLICs 
method and true eUipticity for weight function. Horizontal axis means WSN and vertical 
axis means systematic error. Plots are simulated results of SER of A?^q and line is 2ndSER 
calculated from eq.f lTHj) . 
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Fig. 12. — SER due to RCN for quadrupole moments measured with EGI and E-HOLICs 
method and true eUipticity for weight function. Horizontal axis means WSN and vertical 
axis means systematic error. Cross(square) plots are simulated results of SER of A'Hq(A?^|). 
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Fig. 13. — SER due to RCN for quadrupole moments measured with EGI and E-HOLICs 
method and true eUipticity for weight function. Horizontal axis means WSN and vertical 
axis means systematic error. Square(circle, triangle) plots are simulated results of SER of 
An^iAnl Anf) and line is 2ndSER calculated from eq.dSD, eq.§^, eq.(l83|). 
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From eq.f lMjl A(5(2) is obtained as 



In the situation of EGI, the systematic error is obtained as 



2 



S, 1 - S] 



FigJHlshows the comparison between numerical results of SER and theoretical predation 
of 2ndSER given by eq. (!85l) with <5/ = (0.5,0). The ellipticity is defined by quadrupole 
moments and the averaged effect for quadrupole moments is (eq. (179il and eq. (l80!) ). however 
the ellipticity has non-zero average. It comes from the combination of 1st order effect in 
the complex moments due to I^^^ [6). In this situation, we observe that the ellipticity is 
underestimated by 1.5% in average if WSN = 5. 

We can see the these equations give reasonably good fitting formulas except for the 
images with low WSN. 



4.2. E-HOLICS method with Random Count Noise and without PSF 

Correction 

In the previous section, we derived correction formulas for the systematic error in E- 
HOLICs method in the case of ideal situation, namely we used true ellipticity. In the 
application of E-HOLICs method for realistic situation, we cannot use ellipticity of object 
because of f RCN effect. Thus in reality we define the ellipticities as follows. 



Now we consider RCN effects in the ellipticity for weight function. 
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Fig. 14. — SER due to RCN for observed ellipticity measured with EGI and E-HOLICs 
method with true elhpticity for weight function dw = Sj = (0.5,0). Horizontal axis means 
WSN and vertical axis means SER A(5(2)/5/. Plots are simulated results and line is 2ndSER 
calculated from eq.f l86|) . 



- 30 - 



4-2.1. Centroid Shift Error 

CSE AO is already 1st order of I^'^^{d) and we use square of CSE(i.e. A6'q) in calcula- 
tions of the complex moments. Therefore we may be able to neglect the differences between 
dj and dw in the calculation of CSE because we expect these are higher order effects, Thus 
we obtain 



AO ^ W{Ae,Sw)A0 



and the averaged value as 



2 



rRCN 



Si) 



]9) 



A9 ^W{Ae,Sw)AO = 



A9^^W^{Ae,dw)A9l 



AOl ^ W^{A0,5w)Ael 



w 



WSN^ 1 - 5| 



cr 



W 



WSN^ 1 - 5} 



\C2lS 



C2l 



2Re 



(90) 

C20_C0td)(91) 



^ 5 1 - 2C2°_C0°_ + (CQl\^ 5] j92) 



If the object is EGI, the averaged CSEs are calculated as follows. 



AOl ^ w^{Ae,5w)Ael 



1 



w 



WSN^ 1 - (5| 



AOl ^ W^{Ae,5w)Adl 



w 



WSN^ 1 - 6] 



Si 



Si. 



(93) 

(94) 
(95) 



FigiTS] shows the comparison between numerical results and theoretical prediction given by 
eq.(l93|) and eq.dMD with 6i = (0.5,0). 



4.2.2. Complex Moments 

In this section, we calculate the effects of Ad for the complex moments. 
W {0, Sw) can be expanded up to 2nd order in RCN as 

W{e,Sw) = (6>, 5/ + A5(i) + A5(2)) 

W{e,di){l + ^^ ^1^^ + L«^| (96) 
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Fig. 15. — Plots are A6\Jal^ with elliptical Gaussian image 5i = (0.5,0) and E-HOLICs 
method with observed ellipticity for weight function. Horizontal axis means WSN and verti- 
cal axis means A6'|j/cr^. Square (circle) means A6'q/(T^(A6'|/(T^) measured from simulation 
data and sohd(dashed) line means analytical prediction i.e. eq. fl93|) (eq.f l94|) ). 
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By using this expanded weight, , ^vy) can be expressed by AS as 



^M + A^M(i)('5/) + AHjr,(.)(5/ 



((A5(i) + As^.^ynZl + {^s^,){Sj) + A5(2)(5,)) n^tl) 



{r'^,Z^,Sjl97) 



We postpone to calculate eq.(l97l) in detail until an appropriate expression for AS is available 
later. 



4-2.3. Ellipticity Without PSF Correction 

In this section we calculate the systematic error in measuring the ellipticity, but we 
don't consider PSF smearing. This situation is achieved in space observation such as HST, 
and thus useful. Systematic error in measuring ellipticity of objects with PSF smearing can 
be obtained with more complex calculations and we make a comment in the next section. 

Without PSF smearing, we use the observed ellipticity for the weight function (i.e. 
= Zl, Sw))- By expanding Hi by AS, we obtain 

1 
4 

1 
4 

+-V (A(5(i)AZ«^ + A(5^i) AZ^^) 



)w = nl{r'\zlSw)^nl{r'\zlSi) + AZnl--{AS^,fQl_ + Asi,fQl^^ 



+AH\,) - SjAUl^,^ - - (A(5(2)C0^ + AS^^fOl^) 



_ Ht{r'^ ZlS,) (^^2^^cq4^_ ^ 2|A5(,)rC0t + A5- )COt) 
- {aZ^ - 1 (A5(i)C0^ + ASl.f^lS) 



X 



+ ^ [AS^.^ntir'^ Zl Si) + ASl,{Ht{r'^, Zl Si)) ] ,(98) 



where 



N _ ZZ{I^c\Si) Z^,_,iI^^^,Sj) 
^Zkm= z^,{IoM,Si) Zi{Io^^,Sj) ■ 
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A5(x) are obtained by comparing 1st order as 

^{CAl_AS^i) + COl_ASl)) = AZrI 

A5(i) = 4(cAl*_AZnl-C0l_AZnl*^ 



(100) 



where 



|C42_|2-|C0tP 



coi 



and 2nd order as 



ro^ = 4- 



(101) 



A5(2) = 



An\,^{5i)-SAn\,){Si) 



- \ (A5(2)C0^_ + A5l^fQl_) 
and the average is given by 



^(C4tA5(2) + C0tA(5^2)) 



A^2(2)(5/)-^/A^o(2)(5/)) 



^ 8cr2, (1 - Sj 



WSN^ 



W^) wsF + ^^'^^'-) 

2CAt_ (2 (|C4^_|2 + ICOtp) - |C42_5, - cot 5* 



-C8i2_ (^4C42_C0t + (^4^5/ - COl,S*y^ - C0^_ (^4C42_C0t + (C4:l_di - CO^^*) (l^sj) 

In deriving this expression we used an approximation on the relation of the phase angle(e.g. 
S*nl ^ Snl* see Okura and Futamase 2012). 



In the situation of EGI, eq.( ll03p is reduced to the following expression. 

2 (1 - 6j) 



>w 



AS 



(2) 



^ WSN^ 

2(1-^1) , 
WSN^ ' 



Si 



(104) 
(105) 
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Figure [16] shows the comparison between numerical result of SER and theoretical prediction 
2ndSER given by eq. ffTOSjl . 

By using AS (^2) and eq.( 197|) . the average of A'H^j,^2)i^w) is obtained as 



2a 



w 



1 



8^2 



WSN^ 



alrS 



M 



4WSN2 (1 - 5]\ 

-n^Xl ( 2C42_C0^ (2 - 6]) + (CAl_y 8] + (CO^) ^ (5?^ 



-•Hi^^+l ( 2C4^*_C0t (2 - + (cAl*_) 6 



FigJTT] to FigHni are simulated results of SER of Complex moments. We can see that 
the systematic error of complex moments can be almost estimated by eq.( llU6p expect for 
the objects with low WSN, but there are differences between simulation and theory in the 
estimation of high spin moments "Hi- 



4.2.4- Tests using GREAT08 simulation image 

We test the correction formulas obtained in the previous section using GREAT08 sim- 
ulation data. First, we selected 2 objects from "LowNoise_Known setOOOl.fits" (we call 
Sample A and Sample B), and we compared complex moments and ellipticity between the 
original and noisy objects, where noisy objects are created by adding 10000 types(having 
same RMS) of RCN to the original object. Tabled] and Table [2] are the results of the tests. 
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Fig. 16. — SER due to RCN for observed ellipticity in the situation of EGI and E-HOLICs. 
Horizontal axis means WSN and vertical axis means systematic error A5(2)/<5/. Plots are 
simulated results and line is estimation by eq. dlOSp . 
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Fig. 17. — SER due to RCN for monopole moments measured with EGI and E-HOLICs 
method. Horizontal axis means WSN and vertical axis means systematic error. Plots are 
simulated result of SER of A'Hq and line is estimated value from eq. fllOGp . 
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Fig. 18. — SER due to RCN for quadrupole moments measured with EGI and E-HOLICs 
method. Horizontal axis means WSN and vertical axis means systematic error. Cross(square) 
plots are simulated results of SER A'Hgj^2)(^^2(2)) solid(long dash) line is estimated value 
from eq. fll06p 
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Fig. 19. — SER due to RCN for quadrupole moments measured with EGI and E-HOLICs 
method. Horizontal axis means WSN and vertical axis means systematic error. Square(circle, 
triangle) plots are simulated results of SER of A7{q^2)(^^2(2)' ^"^4(2)) solid(long dash, 
dash) line is estimated value from eq. fllOGp . 
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The 2nd column of the table shows the average of normahzed differences (i.e. the ratio of 
systematic error with and without the corrections) SER and la error and the 3rd column 
shows corrected SER which means systematic error after corrected by 2ndSER and la error 
From these tables, we can see the averaged errors of ellipticity with correction are under 1%. 



Parameter 


SER 


corrected SER 


^0 


0.0386 ±0.00104 


0.00664 ±0.00119 


^0 


0.171 ±0.00426 


0.0118 ±0.00558 


Re[Z|] 


0.233 ±0.00841 


0.0253 ±0.00911 


Im 


0.237 ±0.00832 


0.0189 ±0.0109 




0.491 ±0.0133 


-0.0625 ±0.0235 


Re [Zl] 


0.561 ±0.0172 


-0.0632 ±0.0233 


Im [Zl] 


0.583 ±0.0180 


-0.0876 ±0.0363 


Re [Zj] 


1.24 ±0.224 


-0.286 ±0.556 


Im [Zf] 


0.707 ±0.0235 


-0.0682 ±0.0387 


Re[S] 


-0.0225 ±0.00402 


-0.00636 ±0.00428 


Im [6] 


-0.0199 ±0.00373 


-0.00661 ±0.00400 



Table 1: SER and corrected SER of sample A, average of WSN of noised object is 9.1 and 
ellipticity is (0.39,-0.42) 



5. Conclusion and Future works 

Following the previous work we studied in this paper the systematic error caused by 
signal to noise(SN) ratio of the observed image in our weak lensing analysis called E-HOLICs. 
It has been known that the shear is underestimated when low SN background images are 
used and is overestimated when high SN background images are used in the weak lensimng 
analysis. The latter error was improved in the previous work. The improvement of the former 
error is important because if we have such improvement, we will be able to use many faint 
background sources which improves the statistical accuracy of the weak lensing analysis. 

We have identified the origin of the systematic error as the photon random count noise 
by sky. Although its 1st order effect vanishes by averaging, but 2nd order effects are not 
canceled in measuring the moments and centroid of the images. We investigated this effect 
carefully and obtain the formulas in KSB method and E-HOLICs method to correct the 
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Parameter 


btiix 


corrected SER 


'70 


n noon _i_n nni 1 n 
O.OooO ±O.OOiiO 


n nntxno _i_n nni 1 
0. 00593 ±0.0011o 


'72 


n 1 70 _i_n nn/1/^n 
O.i /^O ±0.004oy 


n nin/i _i_n nn/i^o 
0.0i04 ±0.004M 


T>^ I '721 

Ke [Z2J 


n 001 _i_n nn^/in 
O.Zoi ±0.00(^40 


n ni on _i_n nn^i 
0.0130 ±0.00 / io 


Im [Z2J 


n ootx _i_n ni ntx 
O.zzo ±0.0i0o 


n nntxoi _i_n nnnotx 
0.00531 ±0.009o5 


^0 


U.4yo itU.Uioo 


n nc^fii -un ni on 


Re 


0.551 ±0.0157 


-0.0824 ±0.0152 


Im [Z|] 


0.539 ±0.0191 


-0.107 ±0.0181 


Re [Z^] 


0.648 ±0.0221 


-0.105 ±0.0207 


Im 


0.628 ±0.0229 


-0.142 ±0.0191 


Re [5] 


-0.0254 ±0.00239 


-0.00256 ±0.00252 


Im [5] 


-0.0197 ±0.00650 


-0.00917 ±0.00676 



Table 2: SER and corrected SER of sample B, average of WSN of noised object is 9.1 and 
ellipticity (0.60, -0.23). 

effect in measuring moments and ellipticity. Although general expressions for these formulas 
are complicated, they reduce to relatively simple forms for images with an elliptical Gaussian 
form(EGI). We tested the validity of the correction formula eg. (15 II) and eq. fl63|l for EGIs using 
simulatoin data. Furthermore we applied the general formulas to GREAT08 and confirmed 
that the systematic error reduces to less than 1% in measuring ellipticity for images with 
WSN=9.1 which roughly corresponds to SN=3 object. 

Although the present analysis has not taken into account the Point Spread function(PSF) 
correction which is necessary for the observation from the ground. PSF correction uses com- 
plicated combinations of higher moments and will be very complicated in E-HOLICs ap- 
proach. However the above result is very encouraging and is worthwhile challenging. Finally 
we should point out that the present work will be applicable for the space based observation 
because PSF by instrument is expected to be small for such observation. It will be very 
interesting to confirm this expectation by using data such as COSMOS. 

We thank Satoshi Miyazaki, Takashi hamana, Keiichi Umetsu, Nobuhiro Okabe, Yousuke 
Utsumi and Yuichi Higuchi for useful discussions and comments. This work is partially 
supported by the COE program" Weaving Science Web beyond Particle-matter Hierarchy" 
at Tohoku University and Grant-in -Aid for Scientific Research from JSPS(Nos. 18072001, 
23540282 for TF). 
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A. Appendix 
A.l. GU^w) 

Integration of square Gaussian 
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(Alj 




1 Sw ^ (^w -p^ 


(A2) 
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1 SSwO'w 
~ 4(1 2 


(A6) 
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A. 2. C Coefficients in KSB method with EIHptical Gaussian Image 

In the situation of EGI and KSB method, C coefficients are obtained analyticaUy as 

ci_ = -^(2+^y^2U°'^^o^o) 

C0%^2± = -^mi{r'^,Z,^,0) (A16) 

but foUowing equation cannot be adopt upper general forms 

CX',^ = {X-3-6Sl,)nZ{I'"'',Z,^,Q) (A17) 
where = 0, 2 or 4 and \M\ < \N\. 



A. 3. C Coefficients in E-HOLICs method with true Elhpticity for Weight 

function and Elhptical Gaussian Image 

In the situation of EGI and E-HOLICs method with true elhpticity, C coefficients are 
obtained analytically as 

C^^ = -{l±l){l-6])n%{F'^,Z^,Sj) 

cxZ^ - (x-^-±^^^nZ{r''^,z^,Sr) 

C0^^2± = -^-4^H^Xl{F'\Z,\0) (A18) 

but following equation cannot be adopt upper general forms 

C = -2{l-5])n'o(F'^,Zldr) (A19) 



-43 - 



where iV = 0,2 or 4 and |M| < |A^|. 



nop 

COo = ,^ n lo ^T^ n lo = (A21) 

cot = ,..2 I. '7.02 12 = (A23) 



|C4tP-|C0t| 
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